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Using the spectral function F'(z)/F(z) the RPA correlation energy and other 

properties of a finite system can be written as a contour integral in a compact 

way. This yields a transparent expression and reduces drastically the numerical 



D 
fT ■ efforts for obtaining reliable values. The method applied to pairing vibrations in 



rotating nuclei as an illustrative example. 
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Various problems considered in many-body physics of finite quantum systems can be treated 
in terms of a Hamiltonian which is quadratic in the creation and annihilation operators (see 
for example |jj). The mean field approximation (MF) and the random phase approximation 
(RPA) are the well known examples where the dynamics of complicated system is reducing to 
dynamics described by a quadratic Hamiltonian. These approximations have been proven to 
be quite successful in many applications to mesoscopic systems like nuclei |2|-f|] and metallic 

X 

sj ■ clusters 

. P.. 

One of the long standing problems in the RPA is how to calculate practically the correlation 

energy of the Hamiltonian H = Hq+V, where Hq is a mean field Hamiltonian and V is a two- 
body residual interaction. The solution of this problem could allow also to calculate properly 
other physical quantities which depend on the whole set of the RPA eigen frequencies. 
The MF description is only a first approximation and the RPA correlations can essentially 
improve the description of experimental observables. 

The Hamiltonian, written in terms of the RPA phonon operators 0^(O n ), takes the simple 
form |],| 



H = E MF + J^ ^nO+O n + ^ - E u n Y n+ Y n (1) 

n>0 ^Z 1 n>0 

Here the first term is the MF energy and the second term is the Hamiltonian of the RPA 
(normal) modes. The third term is associated with the contribution of a spurious mode 
arisen due to a broken continuous symmetry (for example translation or rotation). The 
fourth term, where Y n is the vector of the backward-going amplitudes of the RPA solution, 
is a correlation energy e corr caused by the normal modes. The eigen frequencies u u of the 
RPA are derived from the singularities of the RPA response function || 

W(w) = (l-i2°(w)x)- 1 i2°(w), (2) 

where R° is the known unperturbed response function and x is a coupling constant. Ac- 
cordingly, these frequencies are obtained as the zeros of the determinant 

F(u) = det(l-B?(u) X ). (3) 

The total correlation energy E corr , which is essentially the energy gain of the RPA ground 
state relative to the MF ground state, can be expressed as 

J- J corr ^ c-corr = 7T \ / t ^v / t ^n) £-<ex j \f*) 

where E ex is the known exchange energy of the interaction V |1| and E^ are the poles of the 
R°. The eigen frequencies oj v appear in actual situations as rather closely spaced zeros of 
the determinant F{u>), Eq.([3|). Since the number of roots of F(u>) is practically of the order 
10 4 , for example, for heavy nuclei (A fh 150) and none of these roots can be neglected, the 
calculation of the correlation energy E corr is known to be rather difficult even in the case of 
the separable interactions 0]. In addition, in many physical applications it is necessary to 
study the dependence of the correlation energy on the variation of the MF parameters. The 
method, proposed in || for the calculation of the rotational dependence of the correlation 
energy, does not allow treat the contribution of the spurious modes. Furthermore, it should 
be used with a dense grid to approach the necessary convergence in the integration procedure 



to obtain reliable results regarding the contribution of the normal modes. The contour 
integral representation presented in our paper remedies these deficiencies and in addition is 
a transparent method which has further advantages discussed below. 
Let us consider a general eigenvalue equation 

F{z) = (5) 

where F(z) is a supposed to be an analytical function in some region of the complex variable 
z. The complex continuation of the determinant Eq.(^) is an example for an appropriate 
function F(z) studied below in detail. Another example is any finite matrix eigenvalue 
problem F(z) = det (H — z) = written in terms of a finite matrix H, for instance, the 
Hamiltonian of the shell model ||. We define for our purpose the spectral function S(z) as 

using the derivative F'(z) in the z-region where F(z) is analytical. The meaning of S(z) 
becomes obvious for the example F(z) = det (H — z). Assuming H is diagonalizable we may 
express F in terms of the eigenvalues u v as the product F(z) = H u (u u — z) and obtain 

»( \ F ' {z) \~ ! m 

which yields the spectral decomposition of the matrix H. Returning to the general case, it 
is suggestive to apply Cauchy's theorem to the above spectral function S(z). In particular, 
we obtain straightforwardly by forming the following closed contour integral 

V fj, 

where g(z) is an arbitrary complex function which is analytical in the enclosed region. 
Further, the sum v includes all roots uj v of F(z) and the sum \x accounts for all poles p M 
of the derivative F\z) enclosed by the chosen contour. The integers n u and m^ mean the 
multiplicity of the root uj v and the order of the pole p v , respectively. 



The integration formula Eq. (J8|) is now applied to the RPA case. For the sake of the discus- 
sion let us consider a two-body residual interaction which is a sum of separable multipole- 
multipole interactions, i.e. 

v=j:x p q;q p (9) 

p 

Qp=H 1 P ( a f3) a i a ^ + h.c. = Q+ 

a</3 

where a + is a quasiparticle creation operator. This interaction is widely used in nuclear 
physics [0-§] and nowadays it is successfully employed for the description of the plasmon 
frequencies in metallic clusters |nj. The unperturbed matrix response function R° associated 
with the interaction Eq.(^) is defined as the matrix 

fi» = B^ M + ^ M ) do) 

fJL ^P ~ U &fi+U 

where the double index fi = (i,j) is running over all independent quasiparticle pairs i > j = 
l,n and Eft = e« + ej denotes the two-quasiparticle energy. The function F(z) is specified 
as the determinant Eq.(|3[) with matrix elements Eq. (|T0D . The location of its roots u> u is not 
needed but these RPA roots are known to lie on the real axis symmetrically around the origin 
and in addition appear at z — 0, if there are spurious RPA solutions Q. One can easily obtain 
the derivative F'(z) to form S(z). Apparently, the poles of F'(z) coincide with the known 
two-quasiparticle energies E^. For our goal the appropriate contour of the integration Eq.([§D 
goes around all positive roots of F(z) and likewise all two-quasiparticle poles as shown in 
Fig.l. Hence the desired contribution to the RPA correlation energy Eq.(|4|) is obtained as 
the contour integral 

-(yuj u -yE.) = — {dzz^l (ii) 

Usually, the integration will be done numerically. The crucial practical advantage of the 
formula Eq.flTTD is that we are free to choose the rectangular contour (Fig.l) sufficiently 



1 the case of complex RPA solutions is normally excluded but could be similarly treated 



distant from the poles such that the spectral function S(z) = F'(z)/F(z) becomes smooth 
while integrating on this path. 
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FIG. 1. A schematic picture of the integration contour (dashed line) in the complex plane. The 
roots uj v and the poles E^ of the function F{z) are marked with crosses. 

Then, the necessary grid needs not to be dense any more. In practical cases considered so 
far the number of integration points were reduced by a factor 10 2 without loss of a precision 
which reduces the CPU time drastically. The contribution of the possible spurious RPA 
solutions with ujq = is correctly accounted for, since the equivalent sum expression Eq.((4[) 
implies such contributions. Notice that, when integrating in Eq.flll]) separately around such 
a spurious pole at z — 0, there arises no a finite energy contribution due to the z— factor 
in the integrand. The precision of the numerical integration of the correlation energy E corr 
can be conveniently checked by evaluating the analogous contour integral for g(z) = 1 which 
gives according Eq.(J|) an integer number. For the RPA spectral function S(z) the difference 
between the number of RPA roots and the number of quasiparticle states is counted, i.e. 
one expects a zero if there are no spurious uq = solutions. 

We specify discussion to the RPA spectral function S(z) derived from the determinant 
Eqs.(P|),(PH). However, a similar treatment may hold for more general cases. Even, if the 
characteristic function F(z), Eq.(|5|), can be defined only numerically in the vicinity of an 
appropriate integration path one may calculate the spectral function S(z) and apply the 
integration method. Now, we are going to employ the free choice of any physically relevant 
weight function g(z) in the integral formula Eq.(§). In addition, also the integration contour 
can be adapted to the physical problem under study. For g(z) = 1 information on the 
number and distribution of levels can be obtained. First, we define the two-quasiparticle 



spectral function by 



^) = E 7 1 -- (12) 



Let us assume for simplicity, that the derivative F'(z) of the RPA determinant Eqs.(|^), 
(|10|), has simple poles, i.e. m v = 1. Then, in the spectral function 

S(z) = S{z) - S°{z) (13) 

the two-quasiparticle poles are removed. Hence the integral I[g = 1] using the spectral 
function S is counting the number of all RPA states uj v enclosed in the contour. The presence 
of a spurious solution at z = can be checked when calculating T[g = 1] with a narrow 
rectangular contour around zero which cannot be seen directly from the correlation energy 
Eg. ([LTD . The simplest form of a level density distribution can be calculated analogously by 
slices like X[g = 1]/Au; by integrating repeatedly along rectangular stripe contours around 
a given a;— value with an appropriate width Auj. A more refined representation of the level 
distribution follows from Eq.(RI) when inserting the Lorentzian weight function 



A 2 

(u - z) 2 + A 2 



g(z) = L{z) = - - 2 ■ (14) 



which has poles at z = uj ± iA. The RPA spectral function S(z) yields the relation 

D(uJ ' A) ^^ [u - ^y + a 2 = Mm ^ u - * A ) + ^ L ^ ( 15 ) 

where the sum includes all eigenvalues uj v within the chosen contour. Actually it is the 
r.h. side of this relation which can be evaluated numerically and thus it provides us the 
RPA level distribution D(u, A) averaged by the width parameter A. Taking the contour at 
infinity outside all poles (assuming a finite number of roots) the integral becomes zero due 
to the Lorentzian fall off. Hence the integral term in Eq. ([15]) is the Lorentzian " background" 
contribution from the roots not included in the contour. Again, the contour integral involved 
can be calculated conveniently with a smooth integrand as mentioned above. 
To demonstrate the viability and utility of the method, we have applied it to analyze the 
influence of the correlation energy on the description of the g-band of 178 W. For our purpose 



the tilted axis cranking model (TAC) [|TT|] represents the proper Hartree-Fock-Bogoliubov 
(HFB) theory which permits the calculations of the energies and intra band probabilities. 
As a residual interaction we consider the standard pairing field V = —GP + P. It is well 
known that the HFB theory does not provide an adequate description of the transition 
region where the pair field strongly fluctuates. Therefore we use the version of the TAC 
that includes the particle number projection (PNP) and improves the description of the pair 
correlations in high-K bands REg] , Fig.2 shows the results for different Routhian energies 
of the g-band calculated with the TAC alone, with the TAC which incorporates the PNP, 
and the TAC results which include the RPA correlation energy due to pairing vibrations, 
Eq.([|). On the whole region of a rotational frequency the RPA correlation energy lowers 
substantially the TAC Routhian energy and, consequently, improves the description of the 
g-band in comparison with the TAC combined with the PNP. 
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FIG. 2. The g-band of 17S W. Different Routhian energies are presented: for the TAC 
(short-dashed line), for the TAC which includes the PNP (long-dashed line) and for the TAC 
with the RPA correlation energy (solid line). 

Using the same method, we can evaluate the contribution of the RPA correlations to any 
physical quantity in the correlated ground state SQrpa =< Q >rpa — < Q >mf where 
SQrpa = Y,n>o Y,iji ( iij^a*^ji 00- For example, the aligned angular momentum calculated 
in the cranking HFB theory is larger than the experimental alignment. The consideration 
of the additional term arisen due to the RPA correlations leads to better agreement with 



experimental data. Notice that the strength function method [[14] may be considered as an 



alternative approach, however, it would treat the spurious contributions and e corr separately. 
The comprehensive analysis of results obtained with different approaches will be discussed 
in the forthcoming paper |15| . 



In conclusion, the transparent practical method has been developed which drastically re- 
duces the numerical efforts to calculate the integral characteristics renormalized by residual 
interactions in finite quantum systems. In particular, the correlation energy which arises 
in the RPA can be estimated with the necessary accuracy. In the considered example of 
178 W the method demonstrates a remarkable improvement of the description of the g-band 
in comparison with the particle number projection approach. 
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the Heisenberg-Landau program of the JINR. 
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